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ABSTRACT 

We present Karl G. Jansky Very Large Array observations of the CO J = 1 — 0 transition in a sample 
of four z 2 main sequence galaxies. These galaxies are in the blue sequence of star-forming galaxies 
at their redshift, and are part of the IRAM Plateau de Bure HIgh-z Blue Sequence Survey (PHIBSS) 
which imaged them in CO J = 3 — 2. Two galaxies are imaged here at high signal-to-noise, allowing 
determinations of their disk sizes, line profiles, molecular surface densities, and excitation. Using 
these and published measurements, we show that the CO and optical disks have similar sizes in 
main-sequence galaxies, and in the galaxy where we can compare CO J = 1 — 0 and J = 3 — 2 sizes 
we find these are also very similar. Assuming a Galactic CO-to-H 2 conversion, we measure surface 
densities of Emoi ^ 1200 Mq pc“^ in projection and estimate Emoi ^ 500 — 900 Mq pc“^ deprojected. 
Finally, our data yields velocity-integrated Rayleigh-Jeans brightness temperature line ratios rsi that 
are approximately unity. In addition to the similar disk sizes, the very similar line profiles in J = 1 — 0 
and J = 3 — 2 indicate that both transitions sample the same kinematics, implying that their emission 
is coextensive. We conclude that in these two main sequence galaxies there is no evidence for significant 
excitation gradients or a large molecular reservoir that is diffuse or cold and not involved in active 
star-formation. We suggest that rsi in very actively star-forming galaxies is likely an indicator of how 
well mixed the star formation activity and the molecular reservoir are. 
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1. INTRODUCTION 

Galaxy formation and evolution are regulated by the 
interplay between the hierarchical merging of dark mat¬ 
ter halos, the accretion of primordial and recycled gas, 
the transport of gas within galaxy disks, and the subse¬ 
quent gravitational fragmentation, phase transition, and 
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formation of molecular clouds. The star formation that 
is the end-point of this process gives rise to nucleosynthe¬ 
sis, and metal-enriched outffows driven by stellar winds, 
radiation pressure , super novae, and AGN activity (e.g., 
iDave et~aLll2Qlla l ) bl). These mechanisms interact in com¬ 
plex ways, giving rise to the rich phenomenology ob¬ 
served in galaxies today. 

Surveys at z ^ 1 — 2 have shown that most star- 
forming galaxies (SFGs) populate a tight main sequence 
(iNoeske et al.l[2QQ^ iDaddi et al.ll2QQ7 : iRGdighiero et al.l 
suggesting that the evolution of a typical galaxy 
is regulated by secular processes related to gas accre¬ 
tion from the cosmic web, with subsequent star for- 
mation, galactic outffows, and reaccretion c)f gas (e.g., 

Oppenheimer et ^120101: iBonche et al.ll2QlQl:lDave et al.l 

20121: iLillv et al.l 120131: iForbes et al.ll2Q14f ). The gas ac¬ 
cretion is provided by a combination of smooth and 
lumped accretion, the latter likely in the form of mi¬ 
nor mergers (e.g., iKeres et al.l 120051: iDekel et ahl l2009bl: 
iBrooks et al.ll200^ . In this paradigm major mergers play 
a secondary role at forming the majority of galaxies, but 
are still critical for creating the most massive and the 
most actively star-forming galaxies. Until recently, al¬ 
most all of the observations of gas in gal axies at high red- 
shift targeted the rare “behemoths” (see lCarilli fc Wal^ 
l2013[) . typically quasar hosts or merging galaxies bright 
at submi himeter wavelengt hs (or even merging groups of 
galaxies. Ilvison et al.ll2013l ). The exceptions were a few 
intrinsically faint but highly lensed obj ects, where lensing 
complicates the data interpretation (iBaker et al.l l2004 
iGoppin et al.l 120071: iRiechers et al.ll2010[) . The Increased 
sensitivity in the millimeter-wave regime makes it now 
possible to directly observe normal massive star-forming 
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galaxies at high redshift. Indeed, high- 2 ; main sequence 
SFGs are considerably more active in star formation than 
similar mass galaxies at z = 0, likely because they are 
richer in gas than present day blue s equence galaxies, as 
found by molecular gas ob servations (|Tacconi et al.|[2QTQl 
I2Q13I: iDaddi et alJl2QlQallH ). 


Early universe molecular gas surveys usually target the 
mid-J transitions of the CO molecule, which are accessi¬ 
ble for 1 < z < 3 in the transparent 1—3 mm atmospheric 
windows. Additionally, for a fixed gas mass and typical 
conditions these transitions are more luminous and easier 
to detect than the ground level transition. This, however, 
introduces an additional source of uncertaint y when con¬ 
verting their fluxes to H 2 ma sses at high- 2 ; (j Weiss et al.l 
I2QQ7I: ICarilli fc Walter! [2Q13D . Excitation data remain 
scarce for high- 2 ; SEGs. The lowest galaxy luminosities 
currently probed by high redshift CO observations corre¬ 
spond to the high luminosity end of the main sequence, 
and the data suggest their excitation is intermediate be¬ 
tween that observed in high-z Submillimeter Galaxies 
(SMGs) or local Ultr a-Luminous IR Gala xies (ULIRGs), 
and the Milky Way (|Eixsen et al.l 1199^ . Observations 
find CO main-beam velocity-integrated Rayleigh-Jeans 
brightness temperature line ratios r^i = Tmb^v{J = 
3 — 2)/Tmb^v{J = 1 — 0) ^ 0.3 — 0.8 in a handful of 
high-z main sequence SEGs (iDannerbaner et al.l [?009l: 


Arayena et al.l[2QTQt iRiechers et al.l 120101: lAravena et al l 
20141) . The excitation picture, however, is not clear 


cut. Studies of samples of local and high-z luminous 
and ultraluminous IR galaxies also suggests a typical 
rsT ^ 0.6 with very large dispersions (iHarris et al.ll201 


PapadoDoulos et al ] l2011l:llvison et aLll2011l:lGreve et al 


a 




The observed rsi in high-z main sequence SEGs can 
be reproduced by a range of conditions. The typical con¬ 
ditions found in the highly idealized single-component 
models are Tkin ~ 90 K with n ^ 600 cm“^ to Tkin ~ 
10 K with n ~ 2500 c m~^, together with g as filling fac¬ 
tors from 2% to 8% (|Aravena et al.l IMTof ) . This inter¬ 


pretation is by no means unique, however, as realistic 
galaxies contain multiple components, and different tran¬ 
sitions may not be necessarily coextensive — for exam¬ 
ple if dense warm star-forming regions wer e surrounded 
by co oler, less dense envelopes of gas (e.g. IHarris et al.l 
I 2 OIOD . It is also possible to hide a significant amount 
of m ass in a cold component th at would be hard to de¬ 
tect (|Papadoponlos et al.l 12013 ). The inclusion of J > 3 
transitions reveals that, at least in some of the main se¬ 
quence SEGs, the “one component” picture of the ISM 
is too simplistic, although the data support the idea of 
a somewhat lower excitation in main sequence galax¬ 
ies th an in SMGs (|Daddi et al.l |2Q14 iBonrnaud et al.l 
l2Q14al) . In particular, the Tkin and n inferred above 
would produce too little, if any, J > 4 emission in con¬ 
trast with the observations. Better sensitivity observa¬ 
tions of the J = 1 — 0 transition, observations of transi¬ 
tions with J > 3, and higher spatial resolutions are nec¬ 
essary to definitively constrain the molecular excitation 
in these galaxies. 

There are a handful of resolved molec ular measure¬ 
ments of z ^ 1 — 2 main sequence galaxies. iTacconi et al.l 
(IMnh discusses in particular J = 3 — 2 observations of 
one z ^ 1.5 SEG at a resolution of ^ ^ 0'.'65 that reveal 
a collection of clumps with gas masses of ^ 5 x 10^ M 0 , 


sizes of ^ 2—4 kpc, a nd surface densities > 500 Mq pc“^. 
ITacconi et al.l (j2013f ) expands on these results, determin¬ 
ing ^ 20 GO J = 3 — 2 sizes (mostly for z ^ 1 objects), 
reaching the conclusion that t he UV/optica l disk size 
and the GO size are similar. iGenzel et al.l (j2013f ) ex¬ 
amine high-resolution observations of the z ^ 1.5 galaxy 
EGS13011166, a disky system with massive clumps and 
a globally unstable gas disk. lEreundlich et al.l (j2013l ) 
take advantage of position-velocity information to ap¬ 
portion the J = 3 — 2 emission from four z ^ 1.5 SEGs 
in clumps, inferring typical gas surface d ensities Emoi ^ 
400 M e?) pc“^ with a large dispersion. lAravena et al.l 
(|2014f ) presents a few J = 1 — 0 observations of z ^ 1 — 2 
SEGs, including an earlier dataset for the z = 2.2 galaxy 
BX610, one of the sources studied in this paper. In 
the z ^ 1.5 SEG that they observe at high resolution 
they find clumps with masses ^ 4 — 9 x 10^ Mq that 
make up 40% of the emission. Giant clumps are a nat¬ 
ural consequence of the gas richness of thes e high-z sys¬ 
tems, which makes their di s ks unstable (e.g., Dekel et al.l 
l2009al: IGenzel et al.l 120131: IBonrnaud et al.l 2014b[ ). A 
clumpy morphology has also been observed in the rela¬ 
tively rare, very gas-rich and strongly star T erming local 
disks (jGreen et al.ll20l4 lEisher et al.ll2014l) . 

In this work we present new, very sensitive observa¬ 
tions of the CO J = 1 — 0 transition in main sequence 
SEGs at redshift z ^ 2.3, during the peak of cosmic 
star formation. We discuss the observations and data re¬ 
duction in §2. In §3 we present and discuss the results; 
particularly focusing on sizes, surface densities, and ex¬ 
citation. In §4 we summarize them and present the con¬ 
clusions. 


2. OBSERVATIONS AND REDUCTION 

We selected our sources as the brightest CO J = 3 — 2 
emitters at z ^ 2.2 — 2.3 in the IRAM Plateau de 
Bure HIgh-z Blue Sequen ce Survey (PHIBSS) sample 
(|Tacconi et al.l I2Q1Q| . l2Q13f ) . Original optical data for 
these objects are discussed bv lErb et al.l (|2QQ6l ). We ob¬ 
served the CO (J=l — 0) transition at 115.2712 GHz rest 
frequency in Q1700MD94, Q2343BX610, Q1700MD69, 
and Q1700MD174 using the Karl G. Jansky Very Large 
Array (Jansky VLA) in the C- and D-array configu¬ 
rations, with typical angular resolution of ^ 0.7" and 
^ 2" respectively. In order to observe the redshifted CO 
line the observations employ the Ka-band receiver, which 
covers the frequency range 26.5 — 40.0 GHz. 

C-array observations were carried out between 12 Oc¬ 
tober 2010 and 09 January 2011 under project lOB-106. 
These data were obtained during the VLA to Jansky 
VLA upgrade transition period, resulting in a reduced 
bandwidth. Each C-array observation has 0.5 MHz spec¬ 
tral resolution and covers a 128 MHz bandwidth centered 
on the redshifted frequencies of our galaxies (34.55 GHz, 
35.90 GHz, 35.06 GHz, and 34.49 GHz for MD94, BX610, 
MD69, and MD174 respectively). 

The follow up D-array data were taken between 07 
Eebruary 2013 and 07 April 2013 under project 13A-115. 
These data have 1 MHz spectral resolution with a total 
of 4 GHz bandwidth coverage. The time invested is ap¬ 
proximately 8 hours per source in C-array. The D-array 
follow up spent 4.5 hours on MD94 and 8.5 hours on 
BX610. Because only ^ 40% of the visibilities in C-array 
have uv distances < 10^ kA (probing size scales > 2"), 
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and the sources were found to be extended in the orig¬ 
inal observations, the sensitivity added by the D-array 
observations is very significant. 

All data wer e p rocessed using the Jansky VLA re¬ 
duction pipelinJH version 1.2.0 using the Common As¬ 
tronomy Software Applications package (CASA) version 
4.1. The fluxes employed for the flux calibrators were 
applied in the p i peline following the prescribed fits by 
iPerlev fc Bntlerl (j2Q13l ). We visually inspected each ob¬ 
servation after the calibration scripts were completed to 
flag a few remaining bad baselines. The source was then 
split off from each individual observation and combined 
into a single data set for imaging. After imaging with 
50 km s“^ channel widths, we cleaned down to ^ la 
on a 5"x5" box centered on the source. The image 
cubes span a range of ±1000 km s~^ around the op¬ 
tical redshift of the source for MD94 {z = 2.336) and 
BX610 {z = 2.211), while MD69 {z = 2.288) and MD174 
{z = 2.342) were imaged over a ±500 km s~^ span be¬ 
cause of the reduced bandwidth available in the transi¬ 
tion period during which the C-array data were acquired. 
For the calculation of physical parameters we use a stan¬ 
dard cosmol ogy (Hn = 69.6 km s~ ^ Mpc“^, Qm = 0.315, 
^A = 0.685. [Komatsu et al.ll20l3 ). 

BX610: One C-array track was not included in the imag¬ 
ing because of strong pattern noise (maybe due to inter¬ 
ference or correlator problems), but this problem persists 
at a lower level in the remainder of the C-array data for 
this object. The flux and gain calibrators for BX610 
were 3C48 and J2346±0930, respectively. At our observ¬ 
ing frequency 3C48 has a flux density of Sj, = 802 mJy 
which was applied to our data to fix the flux scale. The 
resulting flux of the gain calibrator was 470 mJy in the 
C-array data and 325 mJy in the D-array data. We in¬ 
serted a second calibrator in the D-array observations 
(J2330±110) to check the quality of the calibration in¬ 
cluding phase transfer. The recovered flux of the check 
calibrator is 660 mJy, in excellent agreement with the 
power-law interpolation between the 2 cm and 0.7 cm 
fluxes listed in the VLA calibrator catalog for this ob¬ 
ject of 640 mJy. The synthesized beam sizes range from 
0783 X 0763 (C only, natural weighting) to 1744 x 1721 
(C±D, Briggs weighting with robust parameter set to 0.5, 
henceforth “robust weighting”) and 2702 x 1787 (C±D, 
natural weighting). The la noise per 50 km s“^ channel 
for the robust C±D cube is 52 /iJy beam“^. The high 
resolution J = 3 — 2 data were obtained at PdBI as part 
of PHIBSS, and have a resolution 175 x 077. 

MD94: Along with MD94, flux (3C286) and gain 
(J1645±6330) calibrators were also observed. 3C286 
has a flux density of ?^1800 mJy at 34.55 GHz which 
was used to fix the flux scale of the data. The derived 
flux density of our gain calibrator was 370 mJy in the 
C-array data and 500 mJy in the D-array data. We 
also included a secondary calibrator in the observations 
(J1716±6836). Applying the gain solution to this ob¬ 
ject recovers Sy ~ 430 mJy while the VLA calibrator 
catalog would lead us to expect ^ 670 mJy. However, 
given its listed fluxes this calibrator appears to be a 
highly variable fiat spectrum source, and it is very likely 


that changes are due to intrinsic variability. The synthe¬ 
sized beam sizes range from 0777 x 0772 (C only, natural 
weighting) to 0795 x 0789 (C±D, robust weighting) and 
1770 X 1755 (C±D, natural weighting). The la noise per 
channel for the robust C±D cube is 58 /iJy beam“^. 
MD174 & MD69: These sources were not clearly de¬ 
tected in C-array observations and were not followed with 
D-array. Both galaxies were flux calibrated with observa¬ 
tions of 3C286 and antenna gains were corrected with ob¬ 
servations of J1645±6330. The resulting rms noise per 50 
km s“^ channel is approximately 78 and 68 /iJy beam“^ 
for MD174 and MD69, respectively. The integrated in¬ 
tensity maps of both objects show positive signal at the 
position of the source that, while formally close to 3a 
for MD69, is hard to distinguish from many nearby noise 
peaks. The smoothed spectrum at the position of MD69 
shows what may be a marginally significant line. The 
J = 1 — 0 luminosities of these objects in Tabled] can be 
compared to their respective J = 3 — 2 luminos i ties o f 
1.9 and 1.1 x 10^^ K km s“^ pc^ (|Tacconi et al.ll2QT^ . 
The resulting ratios are consistent with our conclusions 
for MD94 and BX610 in ^3.3[ but do not add significant 
information. 


3. RESULTS AND DISCUSSION 


Both MD94 and BX610 are detected at good signal- 
to-noise, even in the C-only data (Figured])- MD69 and 
MD174, which are weaker sources in the CO 3 — 2 ob¬ 
servations, were not clearly detected and not followed up 
with D array observations, therefore we focus our discus¬ 
sion on the previous two galaxies. Togeth er with the ear¬ 
lier p ublication of another BX610 dataset ([Aravena et al.l 
I2Q14D . these constitute the first CO J = 1 — 0 detections 
of unlensed z > l.b main sequence galaxies, and the first 
observations where measurements of the intrinsic prop¬ 
erties of the cold CO disks of main sequence galaxies free 
of lensing corrections are possible at these redshifts. 

Employing the parametrization of the main sequence of 
star-forming g alaxies derived f o r the AEGIS and COS¬ 
MOS fields bv IWhitaker et aP (j2Q12f ). both MD94 and 


BX610 are very representative of the massive end of the 
main sequence at their z ^ 2.3 redshift. Their specific 
star formation rates are sSER~ 1.8 x 10“^ and 2.1 x 10“^ 
yr“^ respectively, which at their inferred stellar masses 
(M>^ ^ 1.5 X 10^^ and 1.0 x 10^^ Mq respectively) places 
them well within the main sequence 0.3 dex scatter. 
Both are very active star-forming galaxie s, with respec¬ 
tive S ER^ 271±95 and 212±74 Mq yr“^ (|Tacconi et al.l 
l2013[ ). MD94, the most massive of the pair, displays a 
broad Ha line and is classified as an AGN (jErb et al.l 
l2006[ ). BX 610 may contain an AGN, but not a very lumi¬ 
nous one lE d rster Schreiber et al.l l20l4 iNewman et al.l 
l20l4 iGenzel et al .1120 14 ) . It is important to keep in mind 
that these galaxies were selected because they are bright 
GO J = 3 — 2 emitters, which makes it possible to target 
them inGOJ = l — 0, and thus their GO properties may 
not be representative of the population as a whole. Our 
CO measurements are summarized in Table [U and the 
resulting physical param eters in Table [3 toge ther with 
relevant parameters from iTacconi et aH ([201^ . 


The Jansky VLA reduction pipeline is available for down- 3 p The Size of the Moleeular Disks in High-z Main 

load at this url: https://science.nrao.edu/iacihties/vla/data- 
processing/pipeline 
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MD94 is clearly resolved in the observations (c.f., Table 
[T]). The simple Gaussian fit (CASA imfit) yields an 
elongated source in the N-S direction, with deconvolved 
FWHM diameter (after averaging the C-only and C+D 
robust results) of ^ l'/3 x I'/O at high significance (> 5cr). 
With an image scale of ^ 8.4 kpc per arcsecond at this 
redshift, this corresponds to ^ 11 x 8 kpc. The residual 
RMS in the central region after removing the Gaussian fit 
is very similar to that measured in empty regions of the 
original image, showing that a 2D Gaussian is a good 
model for the source at this signal-to-noise. BX610 is 
clearly resolved along its major axis, with a deconvolved 
CO FWHM diameter of 1'.'6 x 0'.'6 corresponding to ^ 
13 X 5 kpc (for this source we adopt the results of the 
C+D robust combination, which despite the larger beam 
is better behaved than the C-only observation). We also 
did UV-plane model fitting using the uvmodelf it task in 
CASA, and the results are essentially identical to those 
from image-plane fitting. 

In z = 0 galaxies there is a good rela tion between 
the s izes of molecular and stellar disks fe.g.. lRegan et al.l 
|2QQ1[) . The CO emission is found in an approximately 
exponential disk with a Rco scale-length that is Rco — 
0.2^25, the isophot al radius of the B-band light at 25^^ 
magnitude (lYonng fc Scovill^ 119911: I Young et al.l 119951: 
iSchrnba et al.ll2Qll[) . Consequently half of the CO emis¬ 
sion in local disks is contained within a radius Ri/ 2 ,co = 
0.34/^25* For optical light it is usually assumed that 
the stellar disk scale length is R^ ^ 0.25^25 (imply¬ 
ing an effective radius containing half the optical light 
of Ri/2,Opt ^ O. 4 R 25 ), so the stellar and molecular disks 
track each other very well in their exponential part. 

How do the CO and optical sizes relate in our high- 
2 : SFGs? For a Gaussian source 50% of the flux is 
contained within the FWHM, therefore the measured 
Ri /2 CO radius for MD94 and BX610 is respectively 
^ 1 / 2 ,CO 4.8 and 4.1 kpc. We can compare directly 
these numbers with near-IR NICMOS observations of 
BX610 (Figure [ 2 ]) that sample the o ptical disk of that 
galaxy (jF orster Schreiber et al.l [2QTT1 ) . Rest-frame 5000 
A observations of BX610 yield an effective radius of 
Ri/ 2,Opt ~ 4.6 kpc, which within the errors is very com¬ 
parable to the J = 1 — 0 Ri/ 2 ,CO we obtain. To put disk 
optical sizes in the 2: ^ 0 context, the typical half-light 
optical radius o f a M:^ ^ 1 0^^ M ^:^ in the SDSS sam¬ 
ple analyzed bv iShen et al.l (j2003f ) is Ri/ 2 ,Opt ^ 4 kpc. 
The existing HST observations for MD94 (Figure [2j) are 
rest-frame 2500 A (jPeter et al.l 120071 ). which combined 
with the AGN presence makes it difficult to obtain a re- 
liable disk size (the 4.8 kpc optical radius reported by 
iTacconi et al.l [20T^ is from Ha). Surface photometry 
shows that 80% of the 2500 A light is in an exponen¬ 
tial component with Ri / 2 ,Opt = Q-37 A 0.08 arcseconds 
(3.1 ± 0.7 kpc) using the lFisher fc Dror^ (j2008[ ) method¬ 
ology. Converting these measurements to an optical size 
is, however, not straightforward eve n in local galaxies 
(which display a range o f scalings, iTavlor et al l 120051: 
iMnhoz-Mateos et al.ll2007[ ) , much less in the poorly char¬ 
acterized high- 2 : population. 

It is interesting to compare our results to some other 
resolv ed measurements of high- 2 : SMGs. iBothwell et al.l 
(l20T0h analyze resolved observations of intermediate ex¬ 
citation CO J = 3 — 2 and J = 4 — 3 in several such 


galaxies at 2 : ^ 1.2 to ^ ^ 2. Their typical CO 50% 
emission radii (computed as the harmonic mean of the 
major and minor axis) are Ri/ 2 ,co ^ 3 kpc, although the 
observations span a range of Ri / 2 ,co ^ f k pc. Similar 
measurements are reported by ITacconi et al.l (j2008[) and 
lEngel et all (j2010[ ). who find typical CO sizes Ri/ 2 ,co ^ 
4 kpc. iSharon et al.l (j2015[) find lensing-corrected source 
sizes of Ri/ 2 ,CO ^ 1-7 kpc in J = 3 — 2 and Ri/ 2 ,co ^ 
2.0 kpc in J = 1 — 0 for the two components of the 
^ ^ 2.7 SMC they study. The Ri/ 2 ,co ^4.1 and 4.8 kpc 
sizes measured here suggest that the CO J = 1 — 0 disks 
of main sequence galaxies tend to be somewhat more ex¬ 
tended than the region producing mid-J CO emission in 
SMGs at comparable redshift. On the other hand, re- 
sol ved CO J = 1 — 0 m easurements of five z ^ 2 SMGs 
by llvison et al.l (j2011f ) find complex morphologies with 
radii of > 6 kpc for four of the objects (the fifth is unre¬ 
solved) and a median FWHM linewidth of 540 km s“^. 
These SMC CO J = 1 — 0 sizes and linewidths are larger 
than those measured for the same objects in the mid-J 
transitions, and also larger than the disk sizes we mea¬ 
sure here. 

The similarity between the molecular and optical 
sizes in main sequence galaxies is also seen in the 
CO J = 3 — 2 measurements reported by ITacconi et al.l 
([ 20 H . mostly obtained for 2: ^ 1 sources (Fig. [3]). 
These measurements suggest that 2: ^ 2.3 main sequence 
galaxy disks follow a scaling between optical light and 
CO J = 1 — 0 similar to that in local disks. The CO is 
not measurably more extended than the stars, as may oc¬ 
cur if a significant amount of pre-enriched newly-arrived 
material is present in the disk outskirts. Such pre¬ 
enriched material is unlikely to arrive through smooth 
accretion; it seems much more likely that occurs through 
mergers, and in particular minor mergers given the place¬ 
ment of these galaxies on the main sequence. The fact 
that the scaling between CO-emitting molecular gas and 
the stars is similar to that observed in local disks sug¬ 
gests that we are not catching these 2: ^ 2 main-sequence 
disks at an out-of-equilibrium stage. Rather, we observe 
them in an equilibrium situation where the established 
stellar population and the enriched CO-bright molecu¬ 
lar gas that fuels the current star formation have similar 
radial distributions. 

3.2. Molecular Masses and Surface Densities 

The excellent signal-to-noise of the observations for 
MD94 and BX610 permits, in principle, a precise de¬ 
termination of their CO J = 1 — 0 luminosity. Unfortu¬ 
nately there is a distressing level of inconsistency between 
the fluxes measured in the C and D configurations, par¬ 
ticularly for MD94, with the former being higher than 
the latter. After detailed inspection of the data and cali¬ 
brations it is not clear what to attribute these differences 
to. In particular, the changes in the flux measured for the 
sources between D and C configurations are not mirrored 
by changes in the flux of their respective gain calibrators 
which could suggest systematic problems with the flux 
calibration (for MD94 the change is in the opposite di¬ 
rection), and are not due to the phase noise as measured 
on the gain calibrator (which shows RMS^ 10 — 20 de¬ 
grees in either configuration). After application of the 
calibrations we also recover the expected flux for the sec- 
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ondary calibrator used in D-array, as noted in §2. The 
flux inconsistency in MD94 is at the ^ 2cr level, which 
suggests that it could be simply due to noise, much higher 
for the C-array data. The fact that the C-array obser¬ 
vations were obtained during a transition period for the 
instrument may also have had an effect on the accuracy 
of those fluxes. Because of these considerations we base 
our discussion on the fluxes from the naturally weighted 
combination of C+D observations, which gives prefer¬ 
ential weighting to the D-array data and has the best 
signal-to-noise. 

To compute t h e mol ecular masses we use Eq. 3 
in iBolatto et al.l (|2Q13f ). which is applicable to CO 
J = 1 — 0, with a Galactic CO-to-H 2 conversion factor. 
This conversion relies fundamentaly on two assumptions: 
1) the molecular gas is primarily self-gravitating, and 2) 
the ratio y/pjT between gas density, p, and tempera- 
ture, T, is roughly s imilar to that in local galaxy disks 
(|Bolatto et al.l 120131 and references therein). We use 
<Tco = 4.36 Mq (K km s~^ pc^)“^ which includes the 
Helium contribution to the mass as is customary (o^co = 
3.2 Mq (K km s“^ pc^)“^ without the Helium contribu¬ 
tion) . A number of arguments point to an approximately 
Galactic value o f the conversion factor for main sequence 
high-^ galaxies ([Tacconi et al.ll2Ql(l iDaddi et al.ll2QlQal: 
iGenzel et al.ll2Q12[ ). including the similarit y between the 
gas m asses inferred from CO and from dust (|Genzel et al.l 
|2Q15[ ). This is supported by detailed numerical model¬ 
ing of high- 2 : galaxies, which suggests that the average 
(Too in these syste ms is not significantly different from 
(Too in local disks (jBournaud et al.ll2Q14a|) . We will see 
in ^3.31 that there are indications that the molecular gas 
temperature is higher in these sources than the average 
in GMCs in the Galaxy. It is important to remark that 
this does not automatically translate into a lower-than- 
Galactic o^co- In a self-gravitating cloud what matters 
is the s/p/T ratio, not just the temperature, as stated 
in assumption number two above. Indeed, there are in¬ 
dications that densities (certainly column densities) are 
much higher in these objects than in Galactic GMCs. 

Under these assumptions, the molecular masses in¬ 
ferred from CO J = 1 — 0 for MD94 and BX610 are 
Mmol ~1.7xl0^^ and 1.1 x 10^^ Mq respectively. These 
masses are sim ilar to the stellar rn asses determined for 
these galaxies (jTacconi et al.l 120131), resulting in molec - 
ular gas fractions of fgas ^ 50%. IBolatto et al.l (|2013f ) 
discuss a surface density correction that is applicable if 
the CO-emitting gas is not self-gravitating but bound to 
the overall potential of the galaxy as is applicable in a 
merger and some galaxy centers (c.f., their Eq. 31). Eor 
our galaxies, application of this correction would amount 
to reducing their molecular mass by a factor of ^ 3. This 
is consistent with what detailed modeling recovers for 
very turbulent starbursting mergers, but it is likely too 
large a correction for the galaxies presented here (by con¬ 
trast modeling of very actively star-forming clumpy main 
seque nce galaxies finds acp within 20% of the Galactic 
value iBournaud et al.ll2Q14aI ) . 

Given the sizes we measure for the molecular emission, 
the molecular masses inferred for MD94 and BX610 im¬ 
ply corresponding average deprojected surface densities 
of Xmoi ^ 900 Mq pc“^ {i ^ 40°) and 400 Mq pc“^ 
{i ^ 65°) assuming the disk size is that of the decon¬ 


volved major axis and the measured elongation is due 
to projection effects (Xmoi ^ 1200 — 1100 Mq pc“^ in 
projection). The surface density found for BX610 is ap- 
pr oximately the typical E moi ^ 400 Mq pc“^ inferred 
by lEreundlich et al.l (j2013f ) using position-velocity infor¬ 
mation in four ^ ^ 1.2 galaxies, as wel l as spatially- 
resolved measurements in 2 : ^ 1.5 galaxies (jTacconi et al.l 
[ 20 T 0 I : IGenzel et al.l [20T1) . MD94, on the other hand, 
exhibits higher deprojected disk surface densities than 
BX610 mostly as a consequence of its smaller inferred 
inclination, that place it at the upper end of the 
lEreundlich et al.l measurements. We note that these “de¬ 
projected” estimates are very tentative and, although 
more physically interesting than projected quantities, 
also considerably more uncertain because of the poor 
knowledge of the geometry. In particular, using the Ha- 
derived inclination of BX610 (see next paragraph) would 
result in a deprojected Emoi ^ 900 Mq pc“^. 

We can also extract size and velocity information from 
the data cubes themselves, and obtain dynamical mass 
estimates. We use GalPak3D, a Bayesian multiparame¬ 
ter Markov Chain Monte Carlo fitter for 3D galaxy data 
that ta kes into account th e effects of instrumental reso¬ 
lution (jBouche et al.ll201^ . to fit the robust C+D cubes 
for MD94 and BX610. The process requires a starting 
set of assumptions (source center, size, inclination, ve¬ 
locity dispersion, etc), that are then used to produce a 
3D model that is compared to the data to compute a 
reduced figure-of-merit which is then minimized. All 
resulting estimates are extremely tentative given the res¬ 
olution and signal-to-noise of the data (Eigure|4]). The 
reduced resulting from the minimization are 1.15 and 
1.05 for MD94 and BX610 respectively, showing that the 
model is a good fit to the data. The minimum, how¬ 
ever, is shallow, showing that the models are not unique. 
BX610 is fit by a (thick) rotating disk with a maximum 
velocity of Vrot ^ 220 km s“^ and intrinsic dispersion 
'^disp ^ 40 km s“^, inclination of i ^ 50°, and radius 
7^1/2,00 ^ 0.4". These parameters are comparable to 
those obtained fro m the analysis of its Ha kinematics 
(jCresci et al.l l2009l Table 2), although both the rota¬ 
tional velocity and radius are lower (the latter is also 
lower than the radius we obtain from the integrated in¬ 
tensity map), and the inclination is higher. The resulting 
dynamical mass computed as M^yn = 2 Ri/2,co'^rot/^ 
using the more robust Ri/ 2 ,co =4.1 kpc estimate from 
Table [ 2 ] is Mdim. ^ IT x 10^^ Mq, lower than found by 
iCresci et al.l (j2009[ ) and too low compared to the inferred 
baryonic mass + M^oi — 2.1 x 10^^ Mq. This is, 
however, entirely attributable to the inclination found 
in the CO fit. If we instead use the Ha- derived in¬ 
clination {i ^ 33°) we correct Vrot to ^ 310 km s“^, 
resulting in M^yn ^ 2.0 x 10^^ Mq in excellent agree¬ 
ment with the estimate of the baryonic mass (we do 
not expect dark matter to make a significant contribu¬ 
tion inside R 1 / 2 ). MD94, on the other hand, is fit as a 
dispersion-dominated system {vdisp ^ 180 km s“^ and 
'^rot ^ 60 km s“^ with i ^ 45°). The corresponding 
dynamical mass is Mdyn ^ 0.9 x 10^^ Mq, too low to ac¬ 
commodate even the stellar mass. Unfortunately there 
are no Ha kinematics to compare with, and the example 
of BX610 highlights the uncertainties of the analysis. 
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An interesting local comparison may be provided by 
the starbursting region of the nearby galaxy NGC253, 
which has 50% of its CO emission within an area of 
150 X 50 pc in projected s ize with a lumino sity Leo ~ 
3.3 X 10^ K km s“^ pc^ (|Lerov et al.l I20l4 Table 5). 
That region of NGC253 has a star formation rate SFR^ 
2 — 3 Mq yr“^, while MD94 and BX610 have inferred 
i ntegrated SFR^ 2 71 and 212 Mq yr“^ respectively 
(jTacconi et al.l l2013f ). About ^ 100 NGC253-hke com¬ 
plexes filling ^ 2% of the area of these main sequence 
disks would be compatible with our CO luminosity and 
surface brightness, while simultaneously fueling the ob¬ 
served star formation activity. This conceptual “picture” 
of an agglomeration of starburst-like regions also fits the 
rsi excitation as we discuss in the next section. Note, 
however, that there is no deep physical significance to 
the inferred beam filling area other than the fact that it 
reproduces the observed integrated surface brightness. 

3.3. The Excitation of CO 

The Cosmic Microwave Background (CMB) has a tem¬ 
perature of Tcmb ~9Katz?^2.3, which sets an 
effective floor for the physical temperature o f the cold 
molecular gas (see e.g.. lda Cnnha et al.ll2Ql^ . In terms 
of the Rayleigh-Jeans brightness temperature (radiation 
temperature) measured by the interferometer, the CMB 
presents a background that lowers the contrast. Thus 
even for optically thick, thermalized emission the ob¬ 
served Rayleigh-Jeans brightness temperature is system¬ 
atically lower than the excitation temperature, Tga^, and 
vanishes w hen = Tomb in the rest frame (see 
Eq. 6 in iBolatto et al.l I2Q13L which is applicable to 
the rest frame). The fact that we see relatively bright 
CO J = 1 — 0 and J = 3 — 2 emission already says that 
Tex > 9 K. 

The effect of the CMB on the ratio of two transi¬ 
tions is small unless the excitation temperature of the 
emission is very low. The Sco^'^ integrated flux ra¬ 
tio between the J = 3 — 2 and J = 1 — 0 transitions is 
10.6 ±1.5 for MD94 and 8.3 ± 1 for BX610. Converting 
to ratios in brightness temperature units requires divid¬ 
ing by the ratio of upper level quantum numbers, J^, 
resulting in rsi « 1.17 ± 0.17 and rsi « 0.92 ±0.11 for 
MD94 and BX610 respectively. These data have much 
better s ignal-to-noise in CO J = 1 — 0 than observations 
used bv lAravena et al.l (j2014l ) to measure rsi = 0.58lo.i3 
in BX610, and the rsi we obtain is at the 1.5cr upper 
boundary of the previous measurement. 

Under the assumption of coextensive emission, these 
integrated numbers suggest that the Tex of the J = 1 — 0 
and J = 3 — 2 transitions are very similar, and both 
transitions arise from optically thick warm gas (that is, 
gas with excitation and kinetic temperature significantly 
higher than the E^ 2 /k ~ 16.6 K energy necessary to col- 
lisionally jump from J = 2 to the J = 3 level, where k is 
Boltzmann’s constant). 

The data strongly suggest that the emission in CO 
J = 1 — 0 and J = 3 — 2 is coextensive in these objects. 
The source sizes are very similar in both CO transitions 
in BX610, where the measured sizes of the CO emis¬ 
sion for the J = 1 — 0 and J = 3 — 2 transitions are the 
same within the uncertainties {R 1/2 = 4.1 ± 1.1 kpc for 
J = 1 — 0, Ri /2 = 3.2 ± 0.6 kpc for J = 3 — 2). More¬ 


over, the spectral shape of the emission is very similar 
for both transitions in MD94 and BX610, within the un¬ 
certainties associated with the noise (Figure [5j). The re¬ 
spective linewidths inCOJ = 1 — 0(COJ = 3 — 2) from 
Single Gaussian fitting are 296 ± 66 (430 ± 61) km s“^ 
and 294 ± 49 (266 ± 35) km s“^ respectively (the dif¬ 
ference in the case of MD94 is driven by the blue side 
of the spectrum, where the fit to J = 3 — 2 incorporates 
as line a wider range of velocities resulting in a shift of 
the centroid). Averaging these results yields a typical 
Av{J = 1 — 0)/Av{J = 3 — 2) = 0.90 ± 0.15 on average 
for these two sources. This suggests that, within our abil¬ 
ity to quantify them with the present data, no dramatic 
excitation gradients are present, as both transitions sam¬ 
ple the same gas kinematics. This is in c ontrast with 
what may be the situation in SMGs, where llvison et all 
(|201lf ) find a measurably larger typical linewidth for the 
COJ=l — Oin comparison to J = 3 — 2, suggesting the 
colder gas is more extended. Note, however, that this is 
a subtle ^ 15% effect in the SMC sample, and so it could 
go undetected in our data. 

These results point to a molecular gas excitation that 
is high and uniform in MD94 and BX610, up to the J = 3 
level. Given the widespread molecular emission, the large 
average gas surface densities, and the inferred SFRs, this 
should not be surprising. The excitation requirements 
for the J = 3 — 2 transition are not particularly strin¬ 
gent: T > Es 2 /k ^ 16.6 K, and an effective critical 
density Uerit ~ 2 x 10^/rco 3-2 ^ few x 10^ cm“^, mak¬ 
ing the rsi ratio a rather blunt indicator of excitation. 
By comparison, the average relation for dust tempera¬ 
ture along th e main sequence pr edicts Tdust ^ 31 K for 
these objects (jCenzel et al.ll2Ql^ Table 5), which can be 
taken as a proxy for the gas temperature in dense en¬ 
vironments and certainly meets the requirement for ex¬ 
citation to J = 3. The measured rsi values show that 
there is not a large reservoir of cold and/or diffuse low 
excitation molecular gas, bright inCOJ=l — 0 but not 
apparent in J = 3 — 2, present in these sources. 

In a study using a l ibrary of galaxy simulations, 
iNaravanan fc Krumholzl (|2Q14h make the case for the 
link between Xsfr and CO line ratios in galaxies. They 
show that in moderately active star-forming galaxies 
(^SFR ^ 0.2 Mq yr“^ kpc“^) the molecular gas has av¬ 
erage temperatures larger than the energy of the J = 3 
level, and given its large typical optical depth it is very 
easy to excite this transition through collisions and ra¬ 
diative trapping. In particular, they parametrize the ex¬ 
pected line ratios as a function of Xsfr- Using their Eq. 
19 and the Esfr ^2 — 3 Mq yr“^ kpc“^ measured for 
our galaxies, we would expect rsi ^ 0.9, very similar to 
what we measure. 

Interestingly, our results are also consistent with the 
rsi measured for CO in the central regions of NGC253 
(jBradford et al.l[2QQ3[) . Another point of comparison is 
provided by the nearby Seyfert NGC1068. There the 
line ratio averages to rsi ^ 1.2 in the starburst ring, 
but with a large range of values correlated with the local 
star formation activity measured in Paschen a\ rsi ^ 
0.7 — 1 in regions wit h low SER, and rsi 2 — 3 in 
regions with high SER (jCarcia-Burillo et al.ll2Ql4) . The 
compact circumnuclear disk region, heavily influenced by 
the AGN, shows rsi ^2.7. 
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Actively star-forming z = 0 galaxies span a range 
of rsi, and although a typical value is rsi ^ 0.66 
there exist several local examples with ^ 0.9 

(Figure [6l iMauersberger et al.l 119991: lYao et all 120031 : 
iPaDadoDoulos et al.ll2011 h There has been a tendency to 
attribute the high rsi measured in some high-^ sources 
to th e presence of a strong AGN (e.g., iRiechers et al.l 
|2on|). Nonetheless, it is clear now that at least some 
high-z s ources display high rsi rat ios with no clear AGN 
activity (|Sharon et al.ll2Ql^ I2015D . That is certainly the 
case also in the NGC253 starburst, where no strong AGN 
is detected, and in the starburst ring of NGG1068 where 
the influence of the AGN is negligible. It is more likely 
that rsi is not so much acting as an indicator of overall 
density and temperature, but rather showing how well 
mixed the star formation activity (which provides the 
bulk of the gas heating) and the molecular reservoir are. 
Galaxies where part of the molecular reservoir is quies¬ 
cent (thus cold and/or diffuse) will show lower rsi, while 
galaxies where most of the molecular gas in experiencing 
star formation will exhibit higher rsi ratios. 

4. SUMMARY AND CONCLUSIONS 

We present and discuss GO J = 1 — 0 observations of 
four main sequence star forming galaxies at 2 : = 2.2 — 2.3 
that are part of the PHIBSS sample. These observa¬ 
tions image and resolve two of these galaxies, MD94 and 
BX610, with very good signal-to-noise. These galaxies 
are representative of the luminous and massive end of 
the main sequence at their redshift. The reader should 
keep in mind that because they were selected on the basis 
of their GO J = 3 — 2 flux, however, it is unclear whether 
they are representative of the CO properties of the wider 
main-sequence population. We use these data to study 
three aspects of these galaxies: 

1. We clearly resolve their molecular disks in CO 
J = 1 — 0, finding half-light CO radii of Ri/ 2,00 ~ 
4.8 and 4.1 kpc respectively. We find that the 
molecular and optical disk of BX610 have half-light 
sizes that agree within the errors, while MD94 lacks 
high resolution rest-frame optical observations to 
compare with . This agrees with the conclusions of 
iTacconi et al.l (|2Q13f) . derived using CO J = 3 — 2 
measurements. It is also similar to the observations 
on local disks where both C O and optical light have 
very similar length scales (| Young fc Scovilld 119911: 
iSchruba et al.|[2QTTI) , suggesting that in these z ^ 2 
main sequence objects we are already observing an 
equilibrium configuration between the established 
stellar population and the CO-bright molecular gas 
fueling star formation. We also observe that the 
CO J = 1 — 0 molecular disk sizes we measure in 
these main-sequence galaxies are more extended 
than those typically obtained for mid-J CO emis¬ 
sion in SMGs at similar redshifts. 

2. We measure molecular masses very similar to 
the stellar masses of these galaxies, and gas sur¬ 
face densities that are 500 — 900 Mq pc“^ depro- 
jected, assuming that the CO-to-H 2 conversion fac¬ 
tor is similar to the Milky Way (which is sup- 
ported in general f or main sequence galaxies, e.g., 
iGenzel et al.ll201§) . The CO surface brightness is 


similar to what would be produced by ^ 100 com¬ 
plexes like that in the center of NGC253 filling 
^ 2% of the area of these disks. This “picture” 
of a collection of NGC253-like starburst complexes 
is also consistent with the observed SFR and rsi 
excitation. 

3. Comparison of our observations with J = 3 — 2 
data yields rsi ratios that are approximately unity. 
This, together with the fact that we measure very 
similar line profiles and disk sizes in both tran¬ 
sitions suggests that these sources do not have a 
large excitation gradient or a significant reservoir 
of cold or low density molecular gas. Because it 
is relatively easy to populate the J = 3 level of 
CO in star-forming gas (where star formation sup¬ 
plies the heating needed to raise the gas tempera¬ 
ture, e.g.. iNaravanan fc Krumholzll2Ql3) . we sug¬ 
gest that rsi may act as an indicator of how well- 
mixed the star formation activity is with the entire 
molecular reservoir. 
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Figure 1. CO J = 1 — 0 integrated intensity maps of our 2 ; ~ 2 main-sequence targets. We show the naturally-weighted C-array 
observations for our four galaxies, and the naturally weighted C+D combination for MD94 and BX610. The corresponding synthesized 
beam is illustrated in the bottom-right corner. The contours start at ±2cr and increase by Icr, with cr listed in Tabled (white contours for 
negative values). The color bar to the right of each panel indicates the intensity scale and the contour values used. 
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Figure 2. Comparison of CO and HST rest-frame optical or UV data. We show the robust Briggs-weighted combination of C+D array 
CO J = 1 — 0 dat a in contours (starting at 2(7, in steps of 2a), overlaid on the HST data. For BX610 the image is HST NICMOS 1.6 jum 
observations from lForster Schreiber et al.l (|2011h . while for MD94 the image is HST ACS 814 nm from IPeter et all (|2007l . kindly provided 
by A. Shapley). The stretch in these images is logarithmic to emphasize the extended emission of the source disks. The corresponding 
synthesized clean beam is illustrated in the bottom right corner. The dashed circle illustrates the i^i/2,co CO J = 1 — 0, our best 
measurement of the intrinsic molecular source size. 



Figure 3. Comparison of optical and CO J = 3 — 2 sizes in main-sequence galaxies (mostly 2 ; ~ 1) reported bv lTacconi et all (|2013I L The 
histogram shows the distribution of the difference between i^i /2 for the optical light, i^i/ 2 ,Opt 5 i^i /2 for the resolved CO measurements 
(mostly from J = 3 — 2), i^i/2,co- The black horizontal bar illustrates the approximate Tier error in the difference (the typical error in 
each size determination is ~ 1.3 kpc). The gray bar shows the J = 1 — 0 measurement presented here for BX610. The distribution of the 
difference of optical and molecular sizes shows that in the vast majority of these galaxies the CO and optical light follow each other, with 
the exception of two galaxies that appear much larger in CO (perhaps undergoing a minor merger). 
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Figure 4. Channel maps showing the J = 1 — 0 emission in BX610 (left) and MD94 (right). We show here the robust combination of 
C+D data. The contours start at 2a and have increments of Icr, with positive values in black and negative in white. Note that the velocity 
interval is 60 km s“^ for BX610 and 100 km s“^ for MD94. The sythesized beams are shown in the bottom right corner. 




Figure 5. Comparison of spatially integrated spectra of MD94 and BX610 in CO J = 1 — 0 and J = 3 — 2. The spectral cubes are 
matched at 2'.'5 resolution for BX610 and 5'' resolution for MD94. The spectral resolution is 50 km s“^ for J = 1 — Q and ~ 28 km s“^ for 
J = 3 — 2. The panels show the spectra for the center position after spatial resolution-matching. The J = 3 — 2 data (in gray) are scaled 
down by a factor of 9 (i.e., J^). 
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Figure 6. Com pilation of rsi values in local samples, measured in nearby galax ies (jManersberger et al.lll99^ L luminous infrared galaxies 
(|Yao et al.ll2003I L and ultra-luminous infrared galaxies (jPapadoDonlos et al.ll 201 ll L All samples have very similar mean values of rsi ~ 0.66. 
The values measured here for MD94 and BX610 (c.f., Table O fall on the high end of the observed range, but are not outliers. 


Table 1 

Photometry 


Dataset 

Position (J2000)*^ 

(7° 

ScoAv d 

^co 

Deconvolved source 

size*^’^ 






Major 

Minor 

P.A. 


RA 

Dec 

(mJy km s“^) 

(mJy km s“^) 

( 10 i°K km s-i pc 2 ) 

(") 

(") 

(°) 

BX610 C 

23*'46“0g?43 

12°49'19'.'3 

14 a 

124 ± 35^ 

3.1 ±0.9 

1.0 ±0.4 

0.4 ±0.2 

27 ±67 

BX610 C±D 

23''46™09!43 

12° 49'1971 

ya 

105 ± 10^ 

2.6 ± 0.2 

1.6 ±0.4 

0.6 ±0.4 

118 ±22 

BX610 CO 3-2 

23*'46™0g?44 

12°49'19'.'3 

40a 

870 ± 63^ 

2.4 ±0.2 

0.8 ± 0.2 

0.7 ±0.2 

34 ± 123 

MD94 C 

00“ 42? 00 

64°11'24'.'3 

14 b 

288 ± 68 ^ 

7.8 ± 1.5 

1.6 ±0.4 

1.0 ±0.3 

6 ±20 

MD94 C±D 

17*'00“42?01 

64°11'24'.'2 

10 ^ 

144 ± 13^ 

3.9 ±0.4 

1.1 ±0.3 

1.0 ±0.3 

114 ±86 

MD94 CO 3-2 

17''00“42!13 

64°11'24'.'6 

170^ 

1520 ± 170^ 

4.6 ±0.5 




MD69 C 

17''00“47?67 

64°09'44'.'3 

10 ^ 

52 ± 18^ 

1.4 ±0.5 

1.3 ±0.6 

0.2 ±0.3 

81 ±28 

MD174 C 



40° 

< 120 ° 

< 3.3° 





^ Integrated over ±150 km s 

^ Integrated between —250 km s“^ and ±500 km s“^ with respect to 2 : = 2.336. 
Integrated over ±300 km s“^, the width of the J = 3 — 2 emission. Limits are 3 ( 7 . 
Deconvolved source FWHM diameters from 2D Gaussian fitting using imfit. 

® In a beam, in the integrated intensity map. 

^ For C±D, reported for robust Briggs weighting 
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Table 2 

Galaxy Parameters 


Galaxy 

Stellar Mass^ 
(10“ Mo) 

SFR'^ 

(M© yr-i) 

Molecular Mass^ 
(10“ M©) 

Diameters^^ 

(kpc) 

■^1/2,CO 
(kpc) 

Ico'^ 

(K km s“^) 

(M© pc-2) 


BX610 

MD94f 

1.0 ±0.3 

1.5 ±0.5 

212 ± 74 
271 ± 95 

1.1 ±0.1 

1.7 ±0.2 

12.2 ±3.4 4.6 ±2.5 
10.8 ±1.2 8.3 ±1.2 

4.1 ± 1.1 
4.8 ±0.4 

248 ± 135 
277 ± 55 

398 ± 146 
928 ± 164 

0.92 ±0.11 
1.17 ±0.17 


^ From ITacconi et al.l (120131) . 

^ Assuming a CO-to-H 2 conversion aco = 4.36 Mq (K km s“^ pc^)“^ which includes the 1.36 Helium correction factor. 

Major and minor diameter enclosing 50% of CO emission computed using imfit. 

^resolved integrated intensity, projected on the sky. 

® Very uncertain deprojected values using size from CO major axis. That has a small effect on MD94, and a much larger effect on BX610. 
Measured “on the sky” values are ~ 1100 — 1200 M© pc“^. 

^ Contains an AGN (|Erb et al.ll2006l) . 
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